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Abstract 

Basic developments in the analytic study of the QCD vacuum 
structure and of the QCD spectrum, including glueballs and hybrids 
are reviewed. 



1 Introduction 

One of the main problems of QCD is the development of quantitative analytic 
nonperturbative methods. Recently a very promising new method has been 
suggested |], (for a review see 0), called the Field Correlator Method 
(FCM) which is fast progressing now. 

FCM has a wide scope and allows to study all nonperturbative problems 
from the QCD vacuum structure [|], || to the details of the QCD spectrum 
||-§|, and high-energy scattering [nj. A practical usefulness of the method 
is connected to the number of gauge-invariant field correlators which should 
be included as the input of the method. Recently the lattice study |TI[ has 



confirmed the phenomenon of the Casimir scaling for the static QQ interac- 
tion, ||, which implies the domimnant role of the quadratic (Gaussian) field 
correlators with accuracy of order of 1%. In this way a very simple picture 
of the QCD vacuum emerges, called the Gaussian Stochastic Approximation 
(GSA), which implies short range correlations between field strength oper- 
ators in the vacuum and the simple spectrum described by the only scale 
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parameter - string tension o (with the a s and correlation length T g defining 
some 10% corrections to the spectrum). 

In this talk we shall review the situation with the QCD vacuum structure 
and the derivation of the QCD spectrum. 

2 The QCD vacuum structure. Stochastic 
vs coherent. 

The basic quantity which defines the vacuum structure in QCD is the field 
correlator (FC) 

£> (n) (zi, -x n ) = (Fmv! (x 1 )$(xi,x 2 )F fi2U2 (x 2 )®(x 2 , x 3 )...F llnl/n (x n )^(x n ,x 1 )), 



The set of FC (H) for n = 2,3, ... gives a detailed characteristic of vacuum 
structure, including field condensates (for coinciding X\ — X2 — ...X n ). On 
general grounds one can distinguish two opposite situations: 1) stochastic 
vacuum 2) coherent vacuum. In the first case FC form a hierarchy with the 
dominant lowest term '(xx,x 2 ) {x\ — x 2 ), while higher FC are fast 

decreasing with n. In the second case all FC are comparable, and expansion 
of physical amplitudes as the series of FC is impractical. This is the case for 
the gas/liquid of classical solutions, e.g. of instantons, magnetic monopoles 
etc. The physical picture behind the situation of nonconverging FC series 
is that of the coherent lump(s), when all points in the lump are strongly 
correlated. 

To understand where belongs the QCD vacuum one can start with the 
Wilson loop in the representation D of the color group SU(3), 



The Stokes theorem and the cluster expansion identity allow to obtain 
the basic equation, which is used in most applications of the FCM (for more 
details see 0) 



(1) 






D (n >(x!, ...x n )da IXlUl (x 1 )...da tlnUn (x n ). (3) 
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Here integration is performed over the minimal surface S m i n inside the 
contour C defined in (|j) and D n is the so-called cumulant or the connected 
correlator, obtained from the FC Eq.(|l|) by subtracting all disconnected av- 
erages. From (g) one easily obtains that the Wilson loop has the area- 
law asymptotics, W{C) ~ exp(— <jS m i n ), for any finite number of terms 
n max ] n < n max in the exponent (g). 

The string tension is expressed through D^ n \ 

° = \j D {2 \x 1 - x 2 )d 2 (x 1 - x 2 ) + 0(flW,n > 4) = a 2 + a 4 + ... (4) 

Eq.(|4]) has several consequences: 1) confinement appears naturally for 
n = 2, i.e. in the GSA 2) the lack of confinement can be due to vanishing of all 
FC, or due to the special cancellation between the cumulants, as it happens 
for the instanton vacuum |4]], 3) for static quarks in the representation D of 
the color group SU(3), the string tension a 2 D ^ is proportional to the quadratic 
Casimir factor ( the Casimir scaling) 

A D) = -^A fUnd ^ C D = \{n 2 + ^ + v 2 + ^ + 3i/). (5) 

However for larger n, n > 4 the Casimir scaling is violated: 

of > = ai Cg J + a 2 (Cf ) 2 + a z cf + ... (6) 

It is remarkable that perturbative interaction of static quarks (r) satisfies 
the Casimir scaling to the order 0(g 6 ) considered so far ||, so the total 
potential V^ D \r) = K,^(r) + cr^ D V + const is also Casimir scaling, if GSA 
works well. 



This picture was tested recently on the lattice |TT] and confirmed the 
Casimir scaling with the accuracy around 1% in the range 0.1 < r < 1.1. fm. 
The full theoretical understanding of this fundamental fact is still lacking, 
both for the perturbative part and for the string tension. On the pedestrian 
level the Casimir scaling and the quadratic (Gaussian) correlator dominance 
implies that the vacuum is highly stochastic and any quasiclassical objects, 
like instantons, are strongly suppressed in the real QCD vacuum. The vac- 
uum consists of small white dipoles of the size T g made of neighboring field 
strength operators. The smallness of T g might be an explanation for the 
Gaussian dominance since higher correlator terms in a are proportional to 
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(FT 2 ) n (T 2 ) -1 , where F is the estimate of the average nonperturbative vac- 
uum field, F ~ 500 (MeV) 2 . 

Lattice calculations of FC have been done repeatedly during last decades, 
using cooling technic [Q and with less accuracy without cooling |L3| . (Re- 
cently another approach based on the so-called gluelump states was exploited 
on the lattice [H| and analytically [15], which has a direct connection to FC). 

The basic result of |12| is that FC consists of perturbative part 0(l/x 4 ) at 
small distances and nonperturbative part 0(exp(— x/T g )) at larger distances 
with T g in the range T g = 0.2 fm (quenched vacuum) and T g = 0.3 fm ( 2 
flavours) . 

Calculations in 
a smaller value 



13 and [14, 15 



16 1 yield 



T w 



take the limit T„ 



as well as sum rule estimates 
0.13 fm to 0.17 fm. This enables us in what follows to 
keeping a = const « 0.18 GeV 2 . 



3 Relativistic dynamics of confinement 

The Gaussian dominance and the small T g limit greatly simplify the rela- 
tivistic dynamics of hadrons. 

As it was shown both for heavy quarks and for light quarks (see || for a 
review) , in the limit of small T g one has local relativistic dynamics which can 
be described by a local Hamiltonian. Before entering into the details of the 
method, it is useful to list a number of problems, present in all dynamical 
calculations of hadrons, which have been easily solved in our Hamiltonian 
method: 

1) the problem of the constituent quark and gluon mass 

2) the problem of Regge slope (a' = ^ for most calculations while in 
nature one has a' = tt—) 

lira ' 

3) the problem of Regge intercept (an arbitrary constant appears in the 
Hamiltonian in most calculations to shift the masses down to the experimen- 
tal values). 

Finally, previous calculations are based on model assumptions and are not 
derived, in contrast to our Hamiltonian directly from the QCD Lagrangian. 
This results in a large number of fitting parameters, whereas in our case we 
have none. 

The basic scheme of the Hamiltonian derivation consists of three steps: 
i) Definition of initial and final states. 
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ii) Exact representation of the gauge-invariant Green's function as a path 
integral over Wilson loops. 

iii) Derivation of the Hamiltonian from the Green's function in the Gaus- 
sian approximation and the small T g limit. 

We start with gauge-invariant asymptotic states, which are built with the 
use of parallel transporters <&(x, y) for nonlocal states. This is done in the 
same way as on the lattice (but without fuzzing links etc.) Thus for the 
meson one has 

*M(x,y) = q{x)r$(x,y)q{y), T = (1, 7^, ...) (g> 1, D v , D^D V ... (7) 

Gluons can be created gauge-covariantly by F^ v and covariant derivatives e.g. 
{D) n F^ v . At this point it is advantageous to separate perturbative (valence) 
and nonperturbative (background) gluonic field B^, |I7| 

A M = 5 M + (8) 

which are integrated independently in the partition function due to the 
'tHooft identity p 



Z=j-J DAe~ s ^ = jf,jDB J Dae~^ B+a \ (9) 

Assigning to homogeneous gauge transformations, one can create valence 
gluonic states with the operators a M , {D) n a^ etc. where D — — igB^. 
In this way glueball fllS] , hybrid ||, [19] and gluelump JT5] states have been 
introduced in good agreement with lattice data. 

As the next step one introduces FFSR for the Green's function ||, for a 
review and earlier refs. see 0]. E.g. for the meson it is 

G r Jx, y) = (G g (, y)TGg(x, y)Y - G q (x, x)YG- q {, y)Y) A (10) 

where FFSR for the (anti) quark Green's function is 

r°° 1 r°° ( dz \ ^ 

G q (x,y) = (m-D)j o ds(Dz) xy e~ K <S> a (x,y), K = - jf dr, (11) 

§a{x,y) = P A exp(ig / A^dz^Pp exp(# / dra^F^) (12) 

Jy JO 

and (Dz) xy implies the path integral from y to x. Here P A , Pf are order- 
ing operators. A similar representation exists for the valence gluon Green's 
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function, with 3 changes; factor (m — D) in front in (|TTD is missing; A^, F^ v 
are in the adjoint representation and ga^F^ in fll2|) should be replaced by 
2gF ap . 

Insertion of ([□]) into (pTOf ) yields the main result to be used below; all 
gluonic field is present in the gauge-invariant form of the closed Wilson loop 
(modulo factor (m — D) and oF insertions which are treated separately). 
Using for the latter the area law, found on the lattice or analytically in the 
Gaussian approximation, one has symbolically (neglecting spin terms) 

G r Jx,y) ~ J dsj Q ds(Dz) xy (Dz) xy e- K - K e-^ {c \ (13) 

Here S m i n (C) is the minimal area inside the contour C made of the variable 
paths C z , C z integrated in (|13D. As the next step one goes over to the Hamil- 
tonian H to escape from the path integrals, using the standard definition 

G r qa (x,y) = (x\e- HT \y). (14) 

Here T is the evolution parameter associated with the chosen hypersurface, 
i.e. for the cm. Hamiltonian it is the cm. time coordinate. To proceed one 



should tackle the proper time s (or r) in (|TT|) , ([131) , an d here enters a new 



fundamental quantity which connects proper and actual (Euclidean) time t, 

dt 

2n(t) = —, t = Zi (r) (15) 



so that the integrals in (0), fljjD can be rewritten as (see for details) 

J ds{D 4 z) xy ... = const J Dfi(t){D 3 z) xy . (16) 



Here enters as an einbein variable P4 , physically the stationary point 
/io of fi(t) will play the role of the constituent mass. Another einbein variable 
v{(3) enters from the representation of the string action in (|TJ) and has the 



physical meaning of the string energy density along the string coordinate 
/3, < (3 < 1. The equal quark mass mi = ni2 = m Hamiltonian was 



obtained in 211 and has the form 



H = P ^±^ + Kr)+ L2/r ' 2 



fi{r) W / u + 2/ 1 (/3-i)M/5)^ 



6 



+ £*V * W 1 ^, (17) 



o z/(/5) Jo 

where p 2 = (pr) 2 /r 2 , and L is the angular momentum, L — (r x p). 

The physical meaning of the terms /i(t) and ^(/?) can be understood when 
one finds their extremal values. E.g. when o = and L = 0, one finds from 



#o = 2^ + m 2 , /i = v / P 2 + ™ 2 ( 18 ) 

so that /xo is the energy of the quark. Similarly in the limiting case L — > oo 
the extremum over z^(/3) yields: 



or 



1 - A = f= , M\L) = 2Kay/L(L + l) (19) 

'1 -% 2 (/5- o) 



so that is the energy density along the string with the /3 playing the role 
of the coordinate along the string. 

A similar form has the two-gluon (glueball) Hamiltonian, which obtaines 
from (|17D putting m = and a — > a ad j = |<r. 

The generalization to more quarks and unequal masses and valence gluons 
is straight-forward, one can find the corresponding examples in H. 



4 Solution of the problems l)-3) 

1) To find the stationary values of [i = /io one can use an approximate relation 
(^|^ =At0 ) pa dE ^) \ fJ - IMj = which has a 5% accuracy 0, and the resulting /i 
which can be considered as a constituent mass, is expressed entirely through 
a. Taking a = 0.18 GeV 2 one obtains for constituent quark mass from 
the lowest meson state nl (meson)=0.35 GeV and for the constituent gluon 
mass from lowest glueball state, /Zq (glueball)= = 0.52 GeV. The proof 
that /io is indeed constituent mass is obtained from the calculation of baryon 
magnetic moments which are expressed entirely through /io in [23 , and agree 
with experiments within 10%. 



2) The Regge-slope problem is actually solved by (|19|), which yield a' = 
at large L due to the third term on the r.h.s. of (0); containing the 
string moment of inertia in the denominator. This is just the term, which 
is absent in the usual relativistic quark model (RQM) predicting a' — 
For not large L the Regge slope (and the physical meson masses) have been 
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calculated in 0, [| in good agreement with experimental meson masses and 



^exp - 



3) the Regge-intercept problem is connected to the fact that actually the 
hadron mass is not given by the eigenvalue Mo of H flT7|), but rather by 
the renormalized value M = Mo — AM. In all RQM calculations AM is 
considered as a phenomenological parameter and is large, AM ~ 0.7 GeV 
for mesons, whereas its physical origin was obscure. 



In |24| it was found that AM appears due to the nonperturbative quark 
mass renormalozation, namely the selfenergy due to the term a^F^ in the 
second order. 

This mechanism gives AM = ^ Y%=i ^Vii where /i, is the constituent 
mass of the i— yh quark, the sum is over all quarks in the hadron and rji(rrii) 
depends only on the current quark mass rrii and ^(0) = 1 for light quarks. 
This expression yields correct values for AM for both mesons and baryons 

0-H; 

It is fascinating that i) /ij, entering the denominator of AM, help to keep 
intact the Regge-slope in M = Mq — AM. ii) the quark mass renormaliza- 
tion 5m 2 q = — is twice as large as in the 1+1 QCD. Thus also the third 
problem - the absolute values of hadron masses - is simply solved within 
the method. The overall agreement of calculated masses of light mesons in 
]7], P , heavy quarkonia in []25| , baryons in || , glueballs in |[J| , hybrids in (19 



with experimental and/or lattice data is better than 10% and is surprising 
since the method exploits only two inputs o and a s (or Aqcd) (i n addition 
to current quark masses renormalized at around 1 GeV). 

Here we have touched only upon two basic topics in the framework of 
the fast developing FCM, more details and applications can be found in the 
review [[| and cited literature. 

The author is grateful to the Organizing Committee of Quarks - 2002 
for kind hospitality and a nice scientific atmosphere at the Conference. The 
financial support of the RFBR grant 00-02-17836 is gratefully acknowledged. 
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